Absence of metallicity in K-doped picene: Importance of electronic correlations 
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Potassium-doped picene (K^picene) has recently been reported to be a superconductor at x = 3 
with critical temperatures up to 18 K. Here we study the electronic structure of K-doped picene fifms 
by photoelectron spectroscopy and ab initio density functional theory combined with dynamical 
mean- field theory (DFT+DMFT). Experimentally we observe that, except for spurious spectral 
weight due to the lack of a homogeneous chemical potential at low K-concentrations (x w 1), the 
spectra always display a finite energy gap. This result is supported by our DFT+DMFT calculations 
which provide clear evidence that K^picene is a Mott insulator for integer doping concentrations 
x = 1, 2, and 3. We discuss various scenarios to understand the discrepancies with previous reports 
of superconductivity and metallic behavior. 

PACS numbers: 71.20.Tx,79.60.Fr,71.15.Mb,71.27.+a 



Since the discovery of superconductivity in layered or- 
ganic charge transfer salts in the late 1970s [l|, a large 
variety of organic superconductors have been reported 
with critical temperatures ranging from a few Kelvins 
to values as high as 38 K in the case of alkali metal- 
doped fullerenes • Whether superconductivity is driven 
by electron-phonon interactions or by electronic correla- 
tions has been a subject of intense debate in the last 
decade, and indeed a significant fraction of these sys- 
tems are presently being discussed as correlated materi- 
als di]. Recently, a pronounced Meissner-Ochsenfeld- 
effect indicating a transition temperature of 18 K was 
observed in polycrystalline powder samples of potassium- 
doped picene [5[ . The parent compound picene is a wide 
band-gap semiconductor (gap ~ 3.3 eV |6j) consisting of 
planar molecules with five aromatic benzene rings (see in- 
set in Fig. [I}. The observation of superconductivity upon 
K intercalation has also been reported in other aromatic 
semiconductors, viz., phenanthrene, coronene and diben- 
zopentacene [7H9J. However, the direct proof of supercon- 
ductivity by reaching zero resistivity is still lacking due 
to the delicate handling of alkali metal-doped organics at 
ambient atmosphere. 

In a simple one-electron picture, K-intercalation is ex- 
pected to induce electron transfer into the lowest unoc- 
cupied molecular orbital (LUMO) of picene, leading to 
a partially filled, i.e., metallic conduction band which 
at low temperatures would host the superconductivity, 
as suggested in various theoretical studies [1, fiol - fisl j . 
However, such calculations include electronic correlations 
only at the level of the exchange-correlation functional 
considered within density functional theory (DFT) and 
neglect electron-phonon interactions. Both effects can 
be quite pronounced in molecular solids [19]. Experi- 
mentally, the picture of the electronic structure has re- 
mained controversial: while one photoemission study re- 



ported the appearance of metallic spectral weig ht at the 
Fermi level for x = 1 [2(| , other authors [2l], [22[ failed to 
observe a metallic Fermi edge for a wide concentration 
range from x — 0.22 to 3.6 and speculated about the 



possible role of correlations and polaronic effects [21 



In view of this unsettled situation we have investigated 
the electronic structure of K x picene in a joint experimen- 
tal and theoretical effort. For this purpose, extensive 
photoemission measurements on pristine and potassium- 
doped picene films have been perfomed, accompanied 
by ab initio density functional theory (DFT) calcula- 
tions in combination with dynamical mean field theory 
(DFT+DMFT). Contrary to the results of Ref.[H both 
experimental and theoretical results confirm the absence 
of metallic behavior for any composition of K^picene with 
integer x = 1, 2, 3. This insulating behavior is attributed 
to strong on-site (i.e., intramolecular) correlations as in- 
ferred from an analysis of the DFT+DMFT calculations. 

Picene films have been grown in ultrahigh vacuum by 
thermal evaporation out of a Knudsen cell on clean na- 
tively oxidized Si substrates. Subsequent doping was 
achieved by evaporation of potassium from commercial 
SAES® dispensers onto the picene films at room tem- 
perature. In situ photoemission spectra were recorded at 
room-temperature with a total energy resolution of 0.1 
and 0.7 eV using He I (21.2 eV) and monochromatized 
Al K Q (1486.6 eV) radiation for valence band and core 
level spectroscopy, respectively. The relative intensity of 
the C Is and K 2p core levels was used to determine the 
doping concentration (see also Supplemental Material). 
All valence band spectra have been corrected for He 1/3 
satellite emission. 

Bulk picene crystallizes in a monoclinic unit cell with 
the molecules arranged in a herringbone structure [23j |. 
As typical for planar molecules on inert substrates, our 
picene films grow in a 3D island mode (Vollmer- Weber 
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type growth), with the long axis standing out from the 
surface in order to minimize film-substrate interaction. 
The c-axis parameter determined from x-ray diffraction 
is f3.51 A, very close to its value in bulk material [El. [23j 
(see Suppl. Mat. for further information). Due to strong 
K-oxidation under ambient atmosphere we have not been 
able to measure reliable structural parameters for the K- 
doped films. 

On the theoretical side the crystal structure of 
K^picene (x = 1,2,3) has been determined by struc- 
tural relaxation within DFT using the generalized gra- 
dient approximation (GGA). The projector augmented 
wave basis was used as implemented in VASP 24427] . 
with additional van der Waals corrections included in 
the relaxation process. The subsequent analysis of the 
electronic structure was performed with the full poten- 
tial local orbital code FPLO [Hj]. 

Figure [lj a) shows the valence band spectrum of a pris- 
tine picene film. It is fully consistent with previous pho- 
toemission studies [U, 0], with the highest occupied 
molecular orbital (HOMO) and following bands clearly 
resolved at -2.2 (H), -2.8 (H-l) and -3.5 eV (H-2) relative 
to the experimental Fermi energy (fj. cxp , determined by 
the Fermi cut-off of a sputtered gold foil) . Figure [TJb) 
displays the evolution of the spectra with K-doping from 
x = to 3 as measured. In agreement with previous work 
the doping results in additional features (Fi, F 2 , F 3 ) in 
the former bandgap and an apparent shift of the overall 
valence band to higher binding energies. This shift is also 
seen in the C Is core level spectrum (not shown here). 

We first discuss the origin of the additional spectral 
features. In analogy to other alkali metal-doped organic 



semiconductors [29|, |30[ F3 is assigned to the relaxed 



HOMO, induced by the structural distortion experienced 
by the doped picene molecules. Its intensity clearly in- 
creases with doping, while the strongly decreasing (but 
still visible) contribution of the H and H-l peaks in Fig.Q] 
reflects the diminishing fraction of uncharged molecules 
upon doping. The weak feature F 2 can be shown to re- 
sult from the oxidation of surface K atoms with oxygen 
present in the residual gas of the ultrahigh vacuum sys- 
tem 31, [iiij]. Therefore, it is a purely extrinsic struc- 
ture not relevant for the following discussions. Finally, 
the most prominent in-gap feature Fi is attributed to a 
(non-rigid-band-like) filling of the former lowest unoccu- 
pied molecular orbitals (LUMOs), as suggested by the 
theoretical calculations discussed below. 

Next we turn to the evolution of the apparent bind- 
ing energies in the valence band spectra. Figure [T] (b) 
shows that for x — 0.8 the new peak Fi sits very close 
or even at the experimental Fermi level, reminiscent of 
the spectrum reported by Okazaki et al. [20] for x = 1, 
while for larger doping concentration it clearly moves 
away from ^ oxp , leaving a (half)gap of ~ 0.5 eV. How- 
ever, the near-Fermi level position of Fi around x — 1 
should not be mistaken for a metallic Fermi edge. First 
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FIG. 1: (a) Valence band spectrum of a pristine picene film. 
H, H-l, and H-2 label the highest occupied molecular orbitals 
of the pristine film, (b) Low-energy part of the valence band 
spectrum showing the evolution with K-doping for various 
near-integer concentrations. The doping-induced features in 
the former bandgap are assigned by Fi, F2, and F3. All en- 
ergies as measured relative to the experimental Fermi energy 
Mcxp- (c) Same spectra, but energetically aligned at the spec- 
tral features of the x = 3 spectrum, for which the existence of 
an equilibrium chemical potential can safely be assumed (see 
text for details). 



of all, its spectral onset is far too broad to account for a 
sharp Fermi edge, even considering thermal and instru- 
mental broadening. Moreover, for smaller concentrations 
between x = and 1 the spectral weight of peak Fi is 
even observed well above the experimental Fermi energy 
(see Fig. S2 in the Suppl. Mat.). Such curious behav- 
ior was observed also in other alkali metal-doped organic 
semiconductor films and was established to result from a 
lack of thermodynamic equilibrium originating from the 
weak intermolecular van der Waals coupling [30( . In such 
a situation, charge fluctuations are strongly suppressed, 
thereby preventing the formation of a well-defined chem- 
ical potential. The ionic charge transfer induced by al- 
kali metal-doping restores the equilibrium. Because in 
our case the doping is achieved by K-evaporation after 
film growth, the doping concentration will decrease from 
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surface to substrate (note that the s-values determined 
from core level spectroscopy are only valid for a thin sur- 
face layer given by the photoemission probing depth of a 
few nm). Only for sufficiently strong overall doping full 
equilibrium within the organic film and between film and 
substrate is achieved. 



In our case we observe saturation of the doping- 
induced energy shifts at x ~ 3, which we identify with 
equilibrium, i.e., the energies of the x — 3 spectrum can 
be safely referred to the experimental Fermi energy. Con- 
sequently, we display in Fig.QJc) all spectra aligned at the 
position of the unrelaxed H and H-l peaks of the x = 3 
film, which incidentally fixes also the doping-induced fea- 
tures at constant energies. This representation nicely 
shows the evolution of in-gap spectral weight upon dop- 
ing, with none of the spectra displaying any metallic 
Fermi edge. The observation that even the x = 3 phase, 
for which superconductivity has been claimed, remains 
gapped is most surpris ing, but in agreement with two 

As one possible explanation po- 
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previous reports 
laronic effects have been suggested [2l| , which may cause 
a pseudogap at the Fermi level as has been observed, 
e.g., in oxides [33j |. However, such a pseudogap typically 
exhibits an exponential decay towards the Fermi energy 
34|, which is not consistent with the present observa- 
tion. Alternatively, the gap may originate from strong 
electronic correlations. In fact, in a Hubbard model-type 
description the intramolecular Coulomb energy U in or- 
ganic solids is often comparable or even l arg er than the 
bandwidth from intermolecular hopping 
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can drive partially filled bands into a Mott insulating 
state. 



In order to further explore the latter scenario, we have 
studied the electronic properties of the DFT-relaxed crys- 
tal structures of K^picene (x = 1, 2, 3) by a combination 
of DFT+DMFT calculations, which include correlation 
effects beyond the GGA approximation in DFT. We start 
first with the DOS for pristine and K x picene obtained 
within GGA (Fig. M (a), (b), (d) and (f)). We observe a 
metallic state for odd x as well as a band gap for x — 2 
as has also been found in previous DFT calculations. 



In DFT+DMFT the kinetic part H K of the Hamilto- 
nian (Eq. [IJ is obtained from the DFT calculations by 
considering the tight-binding approach. Since the DFT 
results indicate some hybridization between the bands 
derived from LUMO and LUMO+1 of the pristine picene 
molecules, we consider for the DMFT calculations a two- 
orbital Hubbard model with a Coulomb interaction U 
and Hund's coupling J z as well as on-site hopping pa- 
rameters i i between molecular orbitals on the same site. 
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FIG. 2: Density of states of pristine and K doped picene. (a) 
DFT DOS for pristine picene with the part of the spectrum 
arising from LUMO and LUMO+1 levels of picene marked by 
gray-shading, (b), (d), (f) DFT DOS of the region close to the 
Fermi level for Kipicene, Epicene and Kapicene, respectively, 
(c), (e) and (g) Combined DFT+DMFT DOS of the same 
compounds. For details see text. 



The Hamiltonian we solve is H = Hk + Hj where 

Hk = - ^ t m C^ m<J C imcr — t± C il(T C l2cr 
{ij)ma i<? 

Hi = U^2 n mir n m ii + ^ [(U - 2J z )n lia n 2 ia+ ( 1 ) 

mi (7 

+ (U - ■iJ z )n 1 i r7 n 2l a-]- 

Bands which are away from the Fermi level are not ex- 
pected to be significantly renormalized by correlations, 
and we save on computational effort by ignoring them 
in the DMFT calculation. As an impurity solver for the 
DMFT equations we employed the weak-coupling con- 
tinuous time Quantum Monte Carlo. Our results are for 
T = 0.025 eV, U = 1.6 eV and J z = (7/8. This choice 
of interaction parameters was guided by recent quantum 
chemistry [361 ] and constrained random phase approxima- 
tion (cRPA) estimates [35J. Note that for Epicene, we 
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find the metal-insulator transition between U = 0.6 eV 
(metallic) and U — 0.8 eV (insulating). The analytic 
continuation of the converged impurity Green's function 
to real frequencies was performed with the maximum en- 
tropy method. The resulting spectral function for these 
two orbitals is shown in Fig. [3J embedded into the DFT- 
DOS of the non-interacting bands. 

Clearly, the DMFT spectral functions show a Mott gap 
not only for doping levels of x — 1 and x = 3, but also for 
x = 2 which in bare DFT is a band insulator with a small 
gap. A closer look at the DOS (Fig. [3] (b)) of K 2 picene 
shows that - contrary to the expectation of a filled picene 
LUMO derived band and an empty picene LUMO+1 de- 
rived band - a Mott gap opens in both orbitals. Inter- 
estingly, energy positions and gap size do not seem to 
depend much on the actual doping level, except for the 
spectral weight transfer from upper to lower Hubbard 
bands, and correspond favorably with the experimental 
spectra, except for some fine structure at higher x not 
resolved in photoemission. The origin of the Mott gap 
for K^picene with x = 1, 2, 3 lies not only in the on-site 
Coulomb interaction U and the Hund's coupling J z but 
also in the interorbital hybridization t±. 

Both experiment and theory thus arrive at the same 
conclusion that K^picene films are (Mott-)insulating for 
all compositions with integer x, including the nominally 
superconducting phase x = 3. This leaves two open ques- 
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FIG. 3: Orbital resolved and total DMFT spectral functions 
of (a) Kipicene, (b) Epicene and (c) Kapicene. 



tions: (1) What is the electronic groundstate for non- 
integer x, and (2) how can the insulating behavior of 
the x — 3 phase be reconciled with the observation of a 
Meissner-Ochsenfeld effect? 

Concerning the first question we note that in a Mott- 
insulator deviation from integer band-filling is normally 
expected to lead to a metallic phase with electron- or 
hole- type carriers [I?} • Note, however, that the doping 
concentration determined in our photoemission experi- 
ments is a macroscopic quantity, averaged over the ex- 
tension of the photon spot (~ 1 mm), and may be com- 
posed of phase-separated integer-filled domains. In fact, 
a recent Raman study 38[ has observed that Kapicene 
decomposes into x = 0, 2 and 3 phases, while a defined 
x = 1 signal could not be detected, indicating that this 
phase is not stable. One can thus understand the pho- 
toemission spectra as a superposition of different integer 
x phases. Because the shape and gapped behavior of 
the corresponding theoretical spectra is very similar, the 
main doping effect will consist in the appearance and 
increasing spectral weight of the LUMO and LUMO+1 
derived lower Hubbard band, as indeed reflected by the 
behavior of feature Fi in our experimental spectra. 

The second question is more puzzling. Our present re- 
sults are in line with the general observation that alkali 
metal-doped organic semiconductors are mostly insulat- 
ing [3(J Hil, with K-doped fullerene (K3C60) being one 
of the rare exceptions [40(. On the other hand, in the 
original paper by Mitsuhashi et al. it has been reported 
that in K^.gpicene the diamagnetic (i.e., superconduct- 
ing) phase is reached from a Pauli paramagnetic, hence 
metallic phase Q . Reconciling our room-temperature re- 
sults with these observations would require an insulator- 
to-metal transition somewhere between 300 and 25 K 
(the highest measured temperature in Ref . [B|) . The Mott- 
Hubbard scenario derived from our DFT+DMFT calcu- 
lations does not provide any mechanism for such a transi- 
tion but would rather stabilize the gap towards low tem- 
peratures. 

Another important aspect is that so far all supercon- 
ducting signatures have been obtained from intercalated 
compact picene powder (HHgf, in contrast to the in vacuo 
grown thin films used here. It is possible that both 
sets of samples differ in subtle structural details account- 
ing for the conflicting findings. On the other hand, the 
DFT+DMFT calculations are based on a DFT-relaxed 
bulk structure which one would expect to be relevant for 
the polycrystalline powder, at least for the interior of the 
powder grains. It is thus tempting to speculate about the 
possibility of the superconducting phase actually being 
located at the grain boundaries, which may have struc- 
tural and/or electronic properties different from the bulk. 
Certainly more clarifying work needs to be done in this 
direction. 

In summary, our combined photoemission and 
DFT+DMFT study finds that Kapicene is a correlated 
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Mott insulator for any integer doping concentration x — 
1, 2 and 3. The insulating behavior of samples with 
macroscopic non-integer filling observed in the experi- 
ment is attributed to microscopic phase separation in 
line with recent Raman results. The apparent spectral 
weight at or even above the Fermi level experimentally 
observed for low doping (x < 1) is attributed to a lack 
of thermodynamical equilibrium, resulting from weak in- 
termolecular van der Waals interactions. The resulting 
discrepancy between the Mott-insulating behavior of par- 
ticularly Kspicene found here and the recent reports of 
superconducting signatures [B} calls for further investiga- 
tions, including the verification of a true non-zero resis- 
tance state in this system. 
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SUPPLEMENTAL MATERIAL 
Structural characterization 

The typical average thickness of the in vacuo grown 
picene films as monitored by a quartz crystal microbal- 
ance (QCM) was 30 nm. The actual film topography as 
measured by atomic force microscopy (AFM) shows that 
the film consists of picene islands with typical lateral di- 
mensions of several 100 nm and heights of up to 80 nm 
(Fig. [Si] (a)). Such 3D island growth mode is expected 
on an inert substrate (here: chemically cleaned and ox- 
idized Si(100)). While the individual islands are crys- 
talline, the overall film appears azimuthally disordered 
with well-developed texture along the c-axis, i.e., perpen- 
dicular to the surface, as inferred from the sharp Bragg 
peaks in the x-ray diffractogram (XRD) in Fig. [Si] (b). 
The resulting c-axis parameter of 13. 51 A confirms that 
the molecules grow in a standing mode. 

Doping concentration 

The atomic K concentration x of the doped picene 
films (Kj; picene) has been determined from the inten- 
sity ratio of the C Is and K 2p photoemission core level 
spectra, taking into account the relative photoionization 
cross sections. The latter is essentially an atomic prop- 
erty and has experimentally been calibrated by inde- 
pendent photoemission measurements of K-intercalatcd 
fullerenes (K3C60), which is an easily reproducible com- 
pound of well-defined stoichiometry. Note that the rel- 
ative cross-section ratio acis/ o~K2p thus determined is 
about 30 % smaller than the value obtained from the 
theoretical cross-sections of Yeh and Lindau [l[. This 
implies that previous photoemission studies which use 
the calculated values overestimate the K concentration 
by the same amount. 



where a homogeneous chemical potential cannot be de- 
fined in the film-substrate system in the early doping 
stages. A possible explanation by a surface photovoltage 
(SPV), often encountered in metal-semiconductor con- 
tacts, can be ruled out, because we do not observe the 
characteristic temperature and intensity dependence of 
the SPV dHl. 

A detailed analysis of identical behavior in an- 
other alkali metal-doped organic semiconductor (para- 
sexiphenyl) has demonstrated that it can be traced back 
to the inability of a purely van der Waals-coupled molecu- 
lar system to establish an equilibrium chemical potential 
throughout the film Q. In such a situation the energy 
alignment refers rather to the vacuum level and thus de- 
pends on work functions. Charge fluctuations induced by 
alkali-metal doping will restore thermodynamical equi- 
librium. Note however, that our potassium deposition 
method leads to a concentration gradient which decays 
from the surface into the film. Low K dosage then leads 
to a doped picene surface layer with little or even 
no K present deeper in the film, maintaining the non- 
equilibrium situation described above. For higher con- 
centrations the film will become doped throughout, hence 
restoring a well-defined chemical potential. The observed 
effect should be less pronounced for thinner films, as is 
indeed the case for the photoemission results of Malms 
et al. 0]. 

We thus take our x = 3 data as an intrinsic equilibrium 
spectrum and have energetically realigned all other spec- 
tra at its spectral features in Fig. IS27 b) , showing gapped 
behavior for all concentrations.. The observation of spec- 
tral weight near the Fermi level reported by Okazaki et 
al. [8( for a nominal doping of x = 1 is attributed to 
the non-equilibrium behavior of peak Fx, which also ex- 
plains why a clear Fermi edge could not be identified in 
that study. 

Structure prediction 



Non-equilibrium effects 

Figure IS2r a) shows a series of valence band spectra as 
measured for a full range of K concentrations from x = 
to 3. The main effect at low doping (x < 1) is a shift of 
the overall spectrum away from the experimental Fermi 
energy /i cxp (which corresponds to the kinetic photoelec- 
tron energy where a metal in electrical contact with the 
sample would display its Fermi edge) and the appear- 
ance of a small peak (Fi) in the former band gap. Closer 
inspection reveals that most of its spectral weight is actu- 
ally located well above the Fermi energy for x — 0.1 and 
only gradually moves below /i cxp with larger doping, until 
the peak finally converges to a constant binding energy 
for x > 1, resulting in a gapped spectrum. This anoma- 
lous behavior clearly signals a non-equilibrium situation 



The relaxation of K^picene candidate structures was 
performed using Fast Inertial Relaxation Engine (FIRE) 
algorithm |9|, emplo yin g Vienna Ab-initio Simulations 



Package (VASP) [lC .lllI w ith the Projector Augmented 
Wave basis (PAW) 



13] in the Generalized Gradient 



Approximation (GGA) as a force and stress tensor calcu- 
lator with plane- wave cutoff of 400 eV and Brillouin zone 
sampled by Monkhorst-Pack /c-point mesh of 5 x 5 x 5 k- 
points. Structures were relaxed until the maximal force 
acting on ions and maximal stress tensor component fell 
bellow 0.01 eV/A and 0.01 GPa, respectively. Fig. [S3] 
shows the structures we obtained. We have compared 
our structures to theoretically obtained lattice parame- 
ters published in Ref. HH Obtained equilibrium volumes 
agree with our results within 2% in all cases, but there 
are some discrepancies in lattice parameters and angles. 
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Without detailed knowledge of the structures in Ref. [14 
it is not possible to say anything further on the degree of 
agreement. 



Details of dynamical mean field theory calculations 

We align Fermi levels of DFT spectra with active bands 
removed and DMFT spectral functions at E = 0. While 
the DFT Fermi level is known for Kipicene and K^picene, 
as they are metallic, and for Kipicene with an error bar 
given by the tiny energy gap of 60 meV, the DMFT chem- 
ical potential follows from the demand that the filling is 
unchanged compared to the DFT calculation. 
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FIG. SI: (a) 5x5/xm 2 AFM image of a nominally 30 nm (QCM) thick picene film showing azimuthally disordered islands, (b) 
X-ray diffraction pattern and rocking curve (inset) of a pristine picene film. The film is highly ordered along the c-axis and 
shows only a small mosaicity. 




FIG. S2: Photoemission spectra of K^picene in a range of ~5eV of the HOMO region for various K-doping concentrations. 
Spectra are shown for x = ... 3. (a) energies as measured, (b) energetically aligned on the spectral features of the x — 3 
compound. 




FIG. S3: Predicted structures of K^picene with x = 1, 2, 3. 



